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Abstract 

Recent molecular dynamic simulations have found chiral single wall carbon nanotubes (SWNTs) 
twist during stretching, which is similar to the motion of a screw. Obviously this phenomenon, as a 
type of curvature-chirality effect, can not be explained by usual isotropic elastic theory of SWNT. 
More interestingly, with larger axial strains (before buckling), the axial strain induced torsion 
(a-SIT) shows asymmetric behaviors for axial tensile and compressing strains, which suggests 
anharmonic elasticity of SWNTs plays an important role in real a-SIT responses. In order to study 
the a-SIT of chiral SWNTs with actual sizes, and avoid possible deviations of computer simulation 
results due to the finite-size effect, we propose a 2D analytical continuum model which can be used 
to describe the the SWNTs of arbitrary chiralities, curvatures, and lengthes, with the concerning of 
anisotropic and anharmonic elasticity of SWNTs. This elastic energy of present model comes from 
the continuum limit of lattice energy based on Second Generation Reactive Empirical Bond Order 
potential (REBO-II), a well-established empirical potential for solid carbons. Our model has no 
adjustable parameters, except for those presented in REBO-II, and all the coefficients in the model 
can be calculated analytically. Using our method, we obtain a-SIT responses of chiral SWNTs 
with arbitrary radius, chiralities and lengthes. Our results are in reasonable agreement with recent 
molecular dynamic simulations. [Liang et. al, Phys. Rev. Lett, 96, 165501 (2006).] Our approach 
can also be used to calculate other curvature-chirality dependent anharmonic mechanic responses 
of SWNTs. 

PACS numbers: 62.25.-g, 46.70.Hg 



2 



The amazing mechanical properties of carbon nanotubes (CNTs), such as high elastic 
modulus, exceptional directional stiffness, and low density, make them idea candidates for 
the applications of nanoelectromechanical systems (NEMS) devises^ 1 ^. Recent studies have 
demonstrated the possibilities of using CNT as actuator 4 , nanotweezers^, and nanorelay^ 1 ^. 
Detailed understanding of mechanical behavior, especially structurally-specific mechanical 
properties of CNT-based NEMS devises is therefore crucial for their potential applications 
in NEMS. 

Unlike isotropic elastic thin shell, due to special geometries of SWNTs, e.g. chiralities, 
there is coupling between axial strain and torsion strain, which is similar to ordinary helical 
spring 10 . More interestingly, recent molecular dynamic simulations found asymmetric behav- 
iors of such coupling in chiral SWNTs 1 *^ and double walled carbon nanotubes (DWNTs) 1 ^, 
namely, asymmetry of a-SIT for tensile and compression strains 1 ^. Later, Upmanyu et a/.'s 
finite element method simulation also obtained asymmetric a-SIT response^ 4 -. Main property 
of asymmetric a-SIT is that a-SIT responses for tension or compression are much different 
at large strain. Torsion angle per unit length increases when strain increases in tension case. 
However, with increasing strain under compression, the torsion angle firstly increases, then 
decreases to zero, and increases again after changing the direction of twistii^ 1 ^. 

A-SIT implies the coupling between axial vibration modes and torsional ones for chiral 
SWNTs, which may play an important role in applications of CNT-NEMS oscillators 1 ^ 1 ^. 
To understand a-SIT response, there are very few studies: Gartstein, et al— used a two- 
dimensional continuum elastic model, predicted linear a-SIT effect for chiral SWNTs with 
small strain, i.e., SWNT twists in opposite directions for tension and compression and ro- 
tation angle varies linearly with strain. Gartstein et al found a-SIT response is chirality 
dependent, it reaches the maximum when the chiral angle is it/ 12. Liang et a/.'s molecular 
simulations 1 ^ extended the study of a-SIT to large strain region (before buckling), obtained 
asymmetric a-SIT. By comparing a-SIT response with changes of geometry of carbon-carbon 
bonds, they found asymmetry a-SIT is relevant to microscopic lattice structure of SWNT. 
Geng et alr^ studied both of torsion induced by axial strain and axial strain induced by tor- 
sion, and showed nonlinear axial stress-strain relation occurring in the same time. Upmanyu 
et ct/'s^ 4 - finite element method simulation also obtained a-SIT. 

All these efforts are valuable in understanding a-SIT. Nevertheless, Gartstein et a/'s theory 
was restrict to linear a-SIT response, while the molecular dynamic or finite element method 
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simulations for a series of SWNTs with some special chiral index were time-consuming, 
a lot of computer resource were needed, which limits their further application to study of 
properties of actual SWNTs. Also, there is a general question for these simulation work: can 
the results of the simulations for small systems be extrapolated to SWNTs at equilibrium 
state with actual sizes? 

In our knowledge, there lacks a easily handled theoretical frame capturing basic physics 
of asymmetric a-SIT which can obtain this response for actual SWNTs at equilibrium states 
with arbitrary radius and chiralities. To fulfill this task, we propose a quasi-analytical 
approach based on continuum elastic theory. In our model, the carbon-carbon interactions 
in SWNT are described by REBO-II potential 17 , which is a classic many-body potential for 
solid carbon and hydrocarbons. The advantages for REBO-II potential are it has analytical 
form of carbon-carbon pair potentials with the bond length and bond angle as variables of 
energy functions, the parameters of REBO-II potential were fitted from a large data sets of 
experiments and ab initio calculations. REBO-II potential can accurately reproduce elastic 
properties of diamond and graphite, In Ref. Illl, molecular dynamic simulation was also based 
on REBO-II potential. 

The carbon-carbon interaction energy near the equilibrium state without deformations 
can be obtained analytically by Taylor expansion with inclusion of the most important cubic 
term, i.e., anharmonic term of bond stretching, 



Here (ij) denotes the nearest neighboring atom pairs, 9^ denotes angle between bonds i — j 
and i — k. Equilibrium state is denoted by "0". Similar series expansion of quadratic terms 
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(1) 



for Brenner potentiate have been reported by Huang et al. 19 . 

The non-crossing second and fourth terms in right hand of Eq. [1] were also presented in 
Lenosky's model 2 ^. From analytical form of REBO-II potential, the derivatives are, 

(%) - 43.67eV-A- 2 , ( f V ) « -5.924eV ■ A -1 , 
\drf J J \dnjd cos 9 ijk J 



d 2 V \ ( d 2 V 

3.187eV, — — — I « -0.367eV, 



d(cos9 ijk ) 2 J ' \d cos 9 ijk d cos 9iji / Q 



and 



d 3 V 
u n / o 



-333.4eV • A" 3 . 



In 2D elastic theory of SWNT, the in-plane deformations of SWNT can be described by 22 




with E\ = En, £2 = £22, ^6 = 2£i2, are the axial, circumferential, and shear strains, 
respectively. After deformation, the bond vector from atom i to its three nearest neighboring 
atoms j, deviates from initial bond vector r$, 

f i:i ss (l + e)^°. 

A SWNT can be viewed as a cylinder with radius R, its surface can be perfectly embedded 
by six-member carbon rings 2 ^. There are three bond curves passing one carbon atoms at 
the surface of SWNT, in the continuum limit, bond vector can be written as 2 ^ 

f°(M) = rg= [l-a 2 /€ 2 (M)/6] a t(M) 

+ [a K(M)/2 + a 2 K s (M)/6] a N(M) (2) 
+ [k(M)t(M)clI/6] a b(M), 

where a = 1A2A is carbon-carbon bond length without strains, M = 1,2,3 denote three 
sp 2 — bonded curves from atom i to atoms j on the surface of SWNT. Vectors t, N and b are 
unit tangential, normal, and binormal vectors of the bond curves from atom % to j, k, r and 
s are the curvature, torsion, and arc parameter of bond curve, respectively, k s = dn/ds.— 
The vectors t(M) = cos 9(M)e x +sm 9(M)e y , b(M) = sin 9 (M) e x — cos 9 (M) e y , where e x and 



e y are the unit axial and circumferential vectors at the z-atom's site on the SWNT surface, 
6(M) is the rotating angle from e x to tangent vector t, which is related to the chiral angle 
9 C M After deforming, bond length = \ fij\, and bond angle between bond vectors and 
are cos 6^ = Uij ■ Uik, with unit vector Uij = fij/rij. Based on these relations, the 2D 
continuum limit of elastic energy per unit area of SWNT in Eq. [TJ which avoids introducing 
ill-defined thickness of SWNTs, can be written as, 

ij i<j<k 

where Cjj and Cijk are in-plane elastic constants, i,j, k = 1,2,6, they have analytical ex- 
pressions, see Appendix. Among them, ciq, harmonic elastic constant for coupling between 
axial strain and torsional twist is proportional to (ao/R) 2 sin(66 l c ), which clearly shows a- 
SIT response is curvature and chirality effect, only occurs in chiral SWNT. Obviously linear 
a-SIT response is distinct at 6 C = n/12 and significant for SWNTs with small diameters, 
which are in accord with previous theoretical and simulation works. For tubes with large 
diameters and small strains, anharmonic elastic energy can be ignored along with c\§ and 
C26 terms, then the isotropic thin shell model for SWNTs is recovered, and the calculated 
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in-plane Young's modulus and Possion's ratio are similar to the results in Ref. 

To study the asymmetric a-SIT, we consider a chiral SWNT with one fixed end, while the 
other end atoms are allowed to relax both radially and tangentially during deformation. The 
axial displacement is fixed for each simulation step, ensuring that only axial stress occurs, 
which is the basic assumption in simulations for a-SIT in SWNTs^ 1 ^ 1 ^. 

The free energy per unit area of SWNT under axial stress is, 

J~ — £elsticity ~ a \ e l- (4) 

Assumption of equilibrium state leads to the following nonlinear equations 

< 5 > 

They give the relation between torsion angle per nm (in unit of degree) <fi = — (I8O/71") x 
Eq/ (i?/lnm) and axial strain e\, which is an asymmetric response. There are two critical 
compressing strains e\ and £**, as shown in Fig{TJ For axial compression, at e*, torsion angle 
reaches its extreme, then SWNT begins to untwist, after totally untwisting at critical strain 
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e**, the tube twists again to the opposite direction, i.e., to the direction as the same as that 
for tension case. 

Another interesting result is nonlinear axial stress-strain relation, the axial secant Young's 
modulus Y s = da j de\ of SWNT is a strict monotonically decreasing function, Y s = Y s0 — 
tei, as shown in Fig. [2J thus SWNTs show strain softening under tension, while strain 
hardening under compression. This phenomenon was also found in recent molecular dynamic 
simulations^. 

We find asymmetric a-SIT and nonlinear axial stress-strain of SWNT are tightly related to 
each other, the nature of which is anharmonicity of atom-atom interaction for SWNTs, such 
as REBO-II potential in Ref. [ll| and present work. This anharmonicity leads to anharmonic 
bond stretching energy in Eq. [T] and cubic terms in Eq. [31 elastic energy. 

To illustrate it, we start from a simplified linear elastic energy per unit area of SWNT, 

F = \ c ^ e \ + c ie £ i £ a + ^c m el - axEx, (6) 

After substituting nonlinear stress-strain relation o\ = Y s0 ei — (t 2 /2) e\ to using equi- 
librium condition dT ' jde\ = 0, the torsion angle, which is proportional to Sq, is a quadratic 
function of axial strain. <ft (e\) curve is a parabola with its symmetric axial located at S\ < 0. 
Thus, present analysis captures main features of asymmetry a-SIT. 

Eq. [3] without cubic terms gives the linear a-SIT response's coefficient 

d4> I _ C12 C26 — C22 Ci6 /_k 

&i 1=0 ~ c 22 c 66 1 j 

with leading term ~ R~ 3 characterizing linear a-SIT response, which is in good agreement 

with Gartstein et aPs theoretical results. Therefore present analysis captures the main 

characters of a-SIT response. 

In our continuum elastic theory, we only get symmetric a-SIT without anharmonic terms, 



however Ref. 



l4 gave asymmetric a-SIT by finite element simulation based on harmonic 



elasticity, although much small (~ 1/1000 of those of molecular simulations). It may be 
due to the elastic energy we used is the continuum limit of lattice energy, which may lose 
some subtle microscopic information. Our continuum elastic theory has some advantages, 
compared to previous simulations, for it is suitable to study SWNTs with actural sizes, and 
all the elastic constants in the theory are obtained analytically. Obviously, our method can 
be extended to calculate other anharmonic properties of SWNTs with arbitrary radius and 
chiralities. 
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In summary, we emphasize the anharmonicity of inter-atoms interactions and curvature- 
chirality induced anisotropic elasticity are both important in a-SIT response, and explain 
the asymmetry a-SIT and nonlinear stress-strain relation all together. We find the unusual 
asymmetric a-SIT effects is the consequence of the curvature-chirality effect and anhar- 
monic elasticity. We give the analytical expressions of anharmonic elastic energy, as well as 
curvature-chirality induced anisotropic elasticity based on REBO-II. The calculated results 
are in reasonable agreement with recent molecular dynamic simulations. Our method can be 
used to analytically calculate anharmonic properties of SWNTs with arbitrary radius and 
chirality. 

We are grateful for helpful discussion with Dr. H. Liang, Prof. Y. Wang and Prof. J. 
Yan. We appreciate Dr. Hangtao Lu for his carefully reading of the manuscripts. 



1 R. H. Baughman, A. A. Zakhidov, and W. A. de Heer, Science, 297, 787 (2002). 

2 H. G. Craighead, Science, 290, 1532 (2000). 

3 S. Sapmaz, Y. M. Blanter, L. Gurevich, and H. S. J. van der Zant, Phys. Rev. B 67, 
235414 (2003). 

4 R. H. Baughman et al, Science 284, 1340 (1999). Philip 

5 P. Kim, and C. M. Lieber, Science 286, 2148 (1999). 

6 J. M. Kinaret, T. Nord, and S. Viefers, Appl. Phys. Lett. 82, 1287 (2003). 

7 S. W. Lee et al, Nano Lett. 4, 2027 (2004). 

8 J. E. Jang et al, Appl. Phys. Lett. 87, 163114 (2005). 

9 J. E. Jang et al, Appl. Phys. Lett. 93, 113105 (2008). 

Y. N. Gartstein, A. A. Zakhidov, and R. H. Baughman, Phys. Rev. B 68, 115415 (2003). 

1 H. Liang, and M. Upmanyu, Phys. Rev. Lett. 96, 165501 (2006). 

2 J. Geng, and T. Chang, Phys. Rev. B 74, 245428 (2006). 

3 H. W. Zhang, L. Wang, J. B. Wang, Z. Q. Zhang, and Y. G. Zheng, Phys. Lett. A, 372, 
3488 (2008). 

4 M. Upmanyu, H.L. Wang, Haiyi Liang, and R. Mahajan, J. R. S. interface, 5, 303 (2008). 

5 V. Sazonova, et al, Nature 431, 284 (2004). 

6 Y. Zhao, C-C. Ma, G. H. Chen, and Q. Jiang, Phys. Rev. Lett. 91, 175504 (2003). 



8 



D. W. Brenner, O. A. Shenderova, J. A. Harrison, S. J. Stuart, B. Ni, and S. B. Sinott, J. Phys. 

Condens. Matter, 14, 783 (2002). 

D. W. Brenner, Phys. Rev. B 42, 9458 (1990). 

Y. Huang, J. Wu, and K. C. Hwang, 74, 245413 (2006). 

T. Lenosky, X. Gonze, and M. Teter, Nature (London), 355, 333 (1992). 

O-Y. Zhong-can, Z-B. Su, and C-L. Wang, Phys. Rev. Lett. 78, 4055 (1997). 

Z-C. Tu, and Zhong-can. Ou-Yang, Phys. Rev. B 65, 233407 (2002). 

Unpublished. 



9 




FIG. 1: Torsion angle-axial strain relations for a series of (8, m) SWNTs, which shows chirality 
dependence of a-SIT response. Only chiral SWNTs have a-SIT response, as shown. 
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FIG. 2: ' Relation between in-plane axial secant Young's modulus and axial strain, for a series of 
(8,m) SWNTs. 
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Appendix 



Elastic constants en, C12, . . . , cqqq presented in Eq. [3] can be described by, 

cc = Mb, 

where, cc is a column vector with the components cc\ to cc\q being the sixteen elastic 
constants of SWCNT, i.e., Cn to respectively, M is a 16 x 6 matrix, b is a column 
vector with components, 




d 2 V \ a , / 8 2 V 
■ —,h 



drijd cos ijk J VLq \d cos 9 ijk d cos 9iji J q fi 



d 3 V 



drf j J Q Q - 

Here, a = 1.42A is carbon-carbon bond length without strains, and = 2.6A 2 is the 
area occupied by one carbon atom at the surface of SWCNTs. 

All 34 non-zero elements of matrix M are analytically written as , 
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A/„,5 = Up} sin(6»), 

Mi « = u + ( is) cos(69) - 

M « = (s) + (is) coa{69h 

= (s) + (il) cos(6e) ' 

Mie,5 = (y^) sin(60). 

Here, a = a /R. 
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